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ABSTRACT 

Aims. The aim of this work is to study the structure of the magnetic field vector in the internetwork and search for the presence of 
small-scale loops. 

Methods. We invert 1.56 fim spectropolarimetric observations of internetwork regions at disc centre by applying the SIR code. This 
allows us to recover the atmospheric parameters that play a role in the formation of these spectral lines. We are mainly interested in 
the structure of the magnetic field vector. 

Results. We find that many opposite polarity elements of the internetwork are connected by short (2 - 6"), low-lying (photospheric) 
loops. These loops connect at least the 10 - 20 % of the internetwork flux visible in our data. Also we have some evidence that points 
towards a dynamic scenario which can be produced by the emergence of internetwork magnetic flux. 
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1. Introduction 

Magnetic fields appear on the solar surface in the form of bipolar 
regions, with the two polarities connected by magnetic field lines 
[ in the shape of loops. If the loops are sufl&ciently large they are 
. traced in Ha (arch filament systems) or EUV or X-ray images 
' (transition region and coronal loops). There is also mounting ev- 
idence for the presence of myriads of small loop-like structures 
in the upper quiet atmosphere (e. g. Feldman et al., 1999) and 
; in the photosphere (iLites et al.Lll996l:[de^Pontieull200 2). and the 
, presence of such loops has been invoked by the complex mix- 
■ ture of magnetic polarities typical of the quiet Sun. Direct mea- 
surements of the full magnetic structure of a loop are extremely 
rare, being basic ally restricted to chro mospheric loops in emerg- 
' ing flux regions (ISolankietal.ll2003h . Here we trace the field of 
low-lying loops and present evidence that they are common in 
the internetwork quiet Sun. 

Our knowledge of the magnetic field in the internetwork 
has evolved rapidly in recent years. The magnetic features 
present in internetwork regions, as seen by the Zeeman efifect, 
are unresolved, occupy a very small portion of the resolu- 
tion element (1 - 2 %) and have field strengths of a few 
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Probably the rest of the resolution element is filled by magnetic 
fields that currently are only detectable by means of the Hanle 
eflTect (LStenflo, 1987) or th e _Zeeman sensitive Mm line at 1.52 
fim dAsensio Ramos et al 1 12007ah . If we assume that such a 
hidden magnetic field is tangled at subresolution scales and that 
it occupies the whole resolution element, the most sophisticated 
Hanle efifect determinations measure a mean magnetic field of 
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100 G ("Truiillo Bueno et al." "2004) while previous estimates 
give a much weaker field (Faurobert et al., 2001). In this type of 
studies, the spatial resolution is given by the slit length, since 
the signal is averaged over the whole slit in order to obtain 
an adequate signal to noise ratio. In additio n, some works 
find empirica l evidence of this tangled nature (StenflgL Il987l: 
Manso Sainz eTall l2004l: lAsen sio Ramos et al, 20071). The 
present paper extends the knowledge of the quiet Sun by reveal- 
ing the connectivity of magnetic features in the internetwork 
and partly its three dimensional structure. 



2. Observations and data reduction 

The observational data consist of two scans of quiet Sun re- 
gions. One of them was located at disc centre and the other 
one at yu = 0.88. They were observed on the 29th and the 30th 
of July 2000, respectively. All four Stokes parameters of the 
pair of Fe i lines at 1 .56 lum were rec orded using the Tenerife 
Infrared Polarimeter dColladosl [199^ installed at the Vacuum 
Tower Telescope at El Teide observatory. In both cases, the inte- 
gration time at each slit position was 1 min which allowed us to 
achieve a noise level in the polarization profiles of 4.1 x 10""^ Ic 
in the observation at disc centre and 2.5 x 10"^ Ic in the one 
at = 0.88. Here, Ic is the continuum intensity. In Fig. [H the 
Ic maps are shown for both scanned areas. The slit was ori- 
ented vertically and scanned horizontally. We used a correlation 
tracker (iBallesteros et al.Lll996l:ISchmidt & KentischeH 1995) to 
correct for image motions due to the seeing. The continuum in- 
tensity contrast of both data sets is 1.4 % and 2.0 %, respectively, 
for the maps at yu = 0.88 and ji = 1 . In both maps, the granulation 
pattern is clear and the spatial resolution was estimated to be of 
the order of T'. This was computed using the Fourier transform 
of the Ic images and choosing the frequency where the power 
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Fig. 1. Continuum maps of the observations taken at disc centre 
(left) and at yu = 0.88 (right). 



was 10 times above noise. In both cases the pixel size, the slit 
width and the scan step were all 0.4''. 

The standard reduction of the data included dark current 
and flatfield corrections. The polarimetric calibration was per- 
formed with the aid of polarization optics of known properties, 
located at an appropriate place in the optical path. Any residual 
crosstalk between the Stokes parameters, because of the uncal- 
ibrated coelostat n iirrors, were rem oved using statistical tech- 
niques described in ' CoUadosI (l2003f) . After that, a statistical pro- 
cedure called Principal Components Analysis (PCA) was ap- 
plied to reduce the noise level of the polarization profiles. This 
consists in creating a basis with the eigenvectors of the correla- 
tion matrix of the observations in which they can be represented. 
Then we can eliminate from the basis those eigenprofiles that 
have no information about the real signals and are dominated by 
noise and reconstruct the observations with this shortened base. 
This is feasible because of the fast decrease of the eigenvalues 
and allows us to increase the noise level of the observations. 

After denoising with PCA, the noise level of Stokes V is 
6.5 X 10"^ Ic and 9.0 x 10"^ Ic for the observations at yu = 1 and 
yu = 0.88, respectively. Stokes Q has a noise level of 1.1 x 10""^ Ic 
and 8.5 x 10"^ Ic and Stokes U 2.0 x lO'^ Ic and 5.0 x 10"^ Ic at 
disc centre and at = 0.88, respectively. 



3. Analysis of the data 

We obtain the magnetic field vector by inverting the four Stokes 
profil es using the SIPQ code ( Ruiz Cobo & del Toro Iniesta, 
Il992h . We perform a two component inversion based on a mag- 
netic atmosphere that occupies a fraction of the resolution ele- 
ment and a field free one filling the rest of the space. Including 
two components in the inversion increases the number of free pa- 
rameters. Then, in orde r to end up with an adequate number of 
degrees of freedom (see lAsensio Ramos et aD. [20Q7bh we chose 
the majority of the variables to be constant with height. The mi- 
croturbulent velocity in both components, the line of sight flow 
velocity of the magnetic component, the magnetic field strength 
and the azimuth of the magnetic field vector are assumed to be 
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constant with height. The macroturbulent velocity is forced to be 
the same in both atmospheres. The macroscopic velocity of the 
non magnetic atmosphere has 3 nodes and the inclination of the 
magnetic field vector has 2 nodes. The temperature stratification 
of each atmosphere was allowed to vary independently with 5 
nodes. 

In Fig. [2] we show the maps of the amplitude of circular po- 
larization for both scans. Superimposed to them, some contours 

of equal amplitude of linear polarization ( + A^) are plot- 
ted (Aq and Au represent the maximum absolute amplitude of 
the Stokes Q and U profiles). In both maps there are numerous 
examples of a linear polarization signal lying in between two 
opposite polarity circular polarization signals. Small low-lying 
magnetic loops are expected to show just such a signature: the 
magnetic field vector points up on one side of the structure, is 
horizontal in the middle and points down on the other side. We 
found 6 clear cases of such structures in the map at disc centre 
and 5 in the map at yu = 0.88. They are enclosed by the white 
rectangles in Fig. |2l In the examples in the yu = 0.88 map all 
the linear polarization is slightly displaced towards one of the 
circular polarization signals. In the disc centre map this is not a 
systematic behaviour. It might be a projection eff'ect or the foot 
points having diff'erent field strengths. 

We discuss now in detail the results from inversions of the 
clearest loop identified by the thick dashed-line rectangle in the 
observation at yu = 0.88 (the others show similar properties and 
geometries). The inclination and azimuth of the magnetic field 
vector recovered from the inversions in the line-of- sight refer- 
ence frame are transformed to the local reference frame (without 
solving the 180° degrees azimuth ambiguity). 

In the left panel of Fig. [3l the horizontal component of the 
magnetic flux vector is presented together with the vertical mag- 
netic flux density for the thick dashed rectangle in Fig. [2l The 
white lines in the left panel of Fig. [3] represent the direction of 
the magnetic field vector. The values of the background flux den- 
sity clearly correspond to internetwork regions in all the pixels 
containing the loop-like structure. The horizontal component of 
the magnetic field vector in the pixels in between the two oppo- 
site polarities is approximately directed along the line that joins 
them. In the central and right panels of Fig.O we plot the mag- 
netic field lines on a vertical cut through the atmosphere along 
the line joining the two polarities. The vertical magnetic flux 
density is indicated by the colour scale and the superimposed 
arrows represent the magnetic field vectors. The azimuth values 
diff'er by 180° in the two panels, both solutions being consistent 
with the observed spectral profiles. Note that the inclination of 
the magnetic field vector is nearly height independent even if we 
allow it to vary. The magnetic field lines have been computed 
as those lines that are parallel to the vector field. In each slit 
position in centre and right panels the geometrical heights have 
been displaced to have the same gas pressure at h=0. The refer- 
ence gas pressure value is the mean value at h=0. The applied 
displacements are small, having a mean value of 20.5 km. The 
Q-loop (centre) or U-loop (right) structure is evident: the mag- 
netic field is almost vertical at the points where the strongest 
Stokes V signals are (at the loop foot points) and is horizontal 
in the transition b e tween o ne polarity and the other. Recently 
iLopez Ariste et al.l (l2006ah have found U-loops in the photo- 
sphere below a filament structure by solving the 1 80° ambiguity 
inherent to the Zeeman eff'ect, suggesting that these structures 
are at least allowed by the physics that drives the solar atmo- 
sphere. In our case, unfortunately, we do not have any observa- 
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Fig. 2. Maps of the amplitude Ay of circular polarization (colour scale) of the observations taken at disc centre (left) and ji = 0.88 

(right). The white contours represent the degree of linear polarization + having the following values: 1.5 x 10"^, 2 x 

10"^, 2.5 X 10"^, 3 X 10"^ Ic for the map at disc centre and additionally at 1 x 10"^ for the map at = 0.88. The white arrow in the 
yu = 0.88 map points towards the position of the disc centre. 



tional constraint to choose between the two different solutions 
showed in Fig.O 

There are many bipolar regions potentially harboring low- 
lying loops in the maps in Fig. [2l but only those boxed in are as 
clear as the example shown in Fig.[3l This means that in all these 
cases the horizontal component of the magnetic vector is mainly 
directed along the line joining the opposite magnetic polarities 
and that the inclination is almost vertical in the foot points of the 
structure and rather horizontal in between. The observed loops 
are rather flat, with a foot point separation of 2300 km and with 
at least some of the field lines at the loop top lying in the height 
range of formation of the 1.56 jim lines (i. e. only a few hundred 
km above the foot points). All cases also share some other prop- 
erties. The magnetic field strength is weak, being below 500 G in 
the majority of the cases. There are only very few points reach- 
ing higher values as 800 - 1000 G. 

The line of sight plasma flows in the magnetic atmosphere 
are correlated with the continuum intensity: bright structures in 
the continuum image are associated with upflows in the magnetic 
component, dark structures with downflows. The foot points of 
the magnetic structures are located equally probably in upflows 
or in downflows as seen in Fig. |4l The more vertical fields are 
located at the footpoints of the structure and we can see that 
they are found either in upflows or downflows. Also the more 
horizontal fields are lying both in upflows and downflows. This 
could be indicative of a dynamic scenario. If the time scale of the 
evolution of such structures is larger than the convective turnover 
time of granulation, one would expect Q-loops to be anchored in 



intergranular lanes (where downflows dominate) and to be lying 
over granules. For U-loops the situation is expected to be re- 
versed, however. Consequently, we may be seeing a mixture of 
Q and U-loops. 

As seen in Fig. [3l the magnetic flux density of the selected 
loop has typical values of the internetwork. This is also true for 
the rest of the loops boxed in Fig.[2l The magnetic flux density in 
the pixels occupied by loops follow the same distribution as the 
pixels in the rest of the internetwork area. Also the magnetic field 
strength distribution is the same as that of the rest of the area cov- 
ered by the typical internetwork. The fact that both the selected 
loops and the rest of the observed area share the same proper- 
ties suggest that the internetwork could harbour many of these 
small low-lying loops, most of them remaining undetectable in 
our observations. 



4. Discussion and conclusions 

We have detected 1 1 examples of low-lying (i. e. photospheric) 
loops in two scans of internetwork regions in the quiet Sun. 
These loops may be related to the horizontal features observed 
by lLites et al.l (Il996h and also the eme rging flux in form of small 
bipoles observed by lde Pontieul (12002 ) using the visible lines at 
630 nm. In our case, we have more sensitive spectropolarimetric 
observations (very magnetically sensitive infrared lines and very 
low noise in both linear and circular polarization) which allow us 
to reconstruct the topology of the magnetic structures. Although 
we could not follow the evolution of such loop-like structures. 
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Fig. 3. Left panel: Map of the vertical magnetic flux density of the loop selected as an example (marked with a dashed square in 
the yu = 0.88 circular polarization map). The white lines represent the horizontal component of the magnetic flux vector. Centre and 
right panels: magnetic field lines computed for the two allowed magnetic field configurations, that diff'er by 180° in azimuth, giving 
rise to an Q- and a U-loop structure, respectively. The inclination and azimuth of the magnetic field vector have been translated to 
the local reference system. 
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Fig. 4. Line-of- sight bulk velocity of the field free component 
plotted versus the inclination of the magnetic field inclination 
with respect to the local vertical. Positive values correspond to 
downflows and negatives to upflows. 



one may speculate that the physical processes that give rise to 
these loops are small scale magnetic flux emergence in the inter- 
network quiet Sun. 

Concerning the magnetic flux carried by the loop-like struc- 
tures that we have identified, we find that at least 6.1% (disc 
centre) and 25.8% (yu = 0.88) of the magnetic flux in the whole 
map is emerging in the form of small scale magnetic loops. On 
average, this means that at least 10-20% of the magnetic flux in 
the solar internetwork is in the form of low-lying loops at any 
given time. This is a lower limit since further such features may 
have escaped detection due to the weakness of the associated 
linear polarization signals. One implication of the fact that these 
loops are so low-lying is that the magnetic flux they are con- 



necting is unlikely to reach the chromosphere and higher layers. 
Unless we are observing a transient phase this implies that pos- 
sibly much of the magnetic flux does not rise above the photo- 
sphere. This point dese rves further atte ntion in view of the sug- 
gestions of Trujillo Bue no et aP (|2004|) that a tangled field could 
provide the clue to understanding how the solar chromosphere 
and corona are heated. 
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